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We demonstrate that modulation of tlie local strength of the cubic self-focusing (SF) nonlinearity 
in the two-dimensional (2D) geometry, in the form of a circle with contrast Ag of the SF coeffi- 
cient relative to the ambient medium with a weaker nonlinearity, stabilizes a family of fundamental 
solitons against the critical collapse. The result is obtained in an analytical form, using the vari- 
ational approximation (VA) and Vakhitov-Kolokolov (VK) stability criterion, and corroborated by 
numerical computations. For the small contrast, the stability interval of the soliton's norm scales as 
AA'' ~ Ag (the replacement of the circle by an annulus leads to a reduction of the stability region by 
perturbations breaking the axial symmetry). To further illustrate this mechanism, we demonstrate, 
in an exact form, the stabilization of ID solitons against the critical collapse under the action of a 
locally enhanced quintic SF nonlinearity. 
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The use of spatially modulated nonlinearities for sup- 
porting solitons in optical waveguides and Bose-Einstein 
condensates (BECs) has drawn much attention in past 
several years [1, 2] . Various patterns of the local strength 
of the nonlinearity can be created in photonic crystals 
or with the help of inhomogeneous distributions 
of resonant dopants in optical waveguides Q- Similar 
structures can be induced in BEC by means of the Fesh- 
bach resonance controlled by nonuniform external fields 
3-[l^. In the latter case, an effective technique for cre- 
ating necessary field patterns may be provided by arrayed 
magnetic films 17]. 

The spatially modulated nonlinearity gives rise to ef- 
fective pseudopotentials [l^ which help to create stable 
spatial solitons [l| in diverse one-dimensional (ID) set- 
tings with the cubic H-Il, quadratic [Ij and 
quintic jisl nonlinearities. Stabilizing multidimensional 
solitons by means of similar techniques is a harder prob- 
lem. In the 2D geometry, smooth landscapes of the self- 
focusing (SF) cubic nonlinearity do not allow one to sta- 
bilize solitons against the corresponding critical ^19'] -[21*1 
collapse [l^. Nevertheless, the stabilization of 2D funda- 
mental solitons in sufficiently broad parametric regions 
is possible in circles [22l . [23| or annuli [13], as well as 
in stripe planforms |24| . with sharp edges. No results 
for stable 2D vortex solitons, or 3D solitons of any kind, 
supported by inhomogeneous SF nonlinearities have been 
reported thus far. On the other hand, the stabilization 
of all kinds of bright solitons and solitary vortices in any 
dimension D is readily provided by the self-defocusing 
cubic nonlinearity, whose local strength is made to grow 
at large distances (r) faster than [2^. A completely 
different mechanism for the stabilization of 2D solitons, 
based on the use of localized gain [l^ , is possible in dis- 
sipative optical media. 



The stabilization mechanism for 2D fundamental soli- 
tons has been demonstrated in areas carrying the SF 
cubic nonlinearity and embedded into self-defocusing or 
linear host media [l^j-jl^l- A challenging question is 
whether the stabilization of solitons against the critical 
collapse may be possible solely due of the contrast be- 
tween stronger SF in a given area, and a weaker self- 
attractive nonlinearity in the ambient medium. The ob- 
jective of this Letter is to demonstrate that such a mech- 
anism works for fundamental solitons trapped in circles 
and annuli. Thus, we demonstrate the stabilization of the 
2D solitons by the enhanced self-focusing: while the well- 
known Townes solitons are unstable against the collapse 
in the uniform SF medium |2H, the creation of a circle 
or annulus with a larger SF coefficient, g, gives rise to a 
family of stable 2D solitons. In particular, in the limit of 
small contrast Ag between the circle or annulus and the 
host medium, the stability interval for the solitons scales 
as A^. To illustrate the genericity of the mechanism, we 
also report exact results for the stabilization of ID fun- 
damental solitons against the respective critical collapse 
induced by the quintic nonlinearity. 

The normalized form of the 2D nonlinear Schrodinger 
(NLS) equation for wave amplitude (x, y, z ) in the 
medium with the modulated nonlinearity is [22| 

= -(l/2)V20-5(r)|0|20, (1) 

where z is the propagation distance, x and y are the 
transverse coordinates, r — \/ x^ + J/^, and 

ofr) = / 1, at p < r < 2, 

■^^ \ 1 - A5, at r < p and r > 2, ^ ' 

p being the inner radius of the annulus (the circle corre- 
sponds to p = 0), while the outer radius is set to be 2 by 
scaling. Inside the annulus, the SF coefficient is normal- 
ized to be 1, while 1 — is the background value. In 
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FIG. 1: (a) A stable localized solution for Ag = 0.5, p = Q 
and — — 0.197, A?^ = 5.73. (b) The propagation constant of 
numerically found families of the trapped fundamental modes 
vs. the norm, at fixed values of the nonlinearity contrast 
between the circle (p = 0) and surrounding background, Ag — 
1.5, 1,0.5, and 0.1. (c) The same as in (b), but as obtained 
from the VA. 



the application to BEC, the evolutional variable in Eq. 
([T|) is time t, instead of z. 

Solutions to Eq. ([l} in the form of axisymmetric local- 
ized modes with propagation constant — /i > are sought 
for a,s (f) — u{r) exp(— z/iz), with real u{r) determined by 
equation 



fiu = -(1/2) (u" + r-'^u') ~ g[r)u 



(3) 



It does not make sense to consider Ag < 0, as in that 
case the annulus is a repulsive structure, that cannot 
trap stable solitons. On the other hand, Ag > 1 implies 
that the host medium is self-defocusing, which, repelling 
the wave field, helps to trap and stabilize modes in the 
SF annulus [1^. In the case of Ag > 1, the field is 
strongly confined to the annulus. Indeed, in a boundary 
layer adjacent to r = 2, the solution to Eq. ([3]) interpo- 
lating between the inner one (at r < 2), which is char- 
acterized by the boundary value of the derivative, uJ,^2j 
and the outer solution, which quickly decays at r > 2, is 

'{Agf* + v/A^ 



^^- 



with 



M (r = 2) w {Ag)~''''^ •\/~'"r=2- This means that, in the 
case of Ag — oo, one should look for the inner solution 
satisfying the boundary condition u{r = 2) = 0, with 
arbitrary The same pertains to the boundary at 

r = p. 

Numerical solutions of Eq. ([3]) were found by means 
of the shooting method. Figure [IJa) displays a typical 
profile of stable fundamental (single-peak) modes. Fam- 
ilies of the solutions are presented in Fig. [Ijb) by curves 
for the total power (norm), iV = 27r v?{r)rdr, ver- 
sus /i. All the families contain portions that meet the 
Vakhitov-Kolokolov (VK) stability criterion, d^/dN < 
[Ullil,!!^. Direct simulations of the full 2D equation (P) 
(without imposing the axial symmetry) confirm that all 
the solutions satisfying the VK criterion are stable, while 
those corresponding to dfi/dN > suffer the collapse or 
decay, as expected. 

A natural tool for the consideration of the present 
setting is provided by the variational approximation 
(VA) (2^, [2^. To this end, we adopt the Gaus- 
sian ansatz, u = ^ exp{— r^/(2W^^)}, with amplitude 



A, width W, and norm N = nA^W^. The substi- 
tution of this into the underlying Lagrangian, L = 

27T 2/iM^ — (u')^ + g{r)u'^ rdr, yields 



N 



[1- Ag(l-i?2+i?i)], (4) 



where Ei = exp(-8/W^2) ^nd E2 = ey.^{-2p^ /W"^). 
Then, W and ^ are determined by the variational equa- 
tions, dL/dW = dL/dN = 0. FigurefHc) shows that the 
VA produces curves fi{N) which are very close to their 
numerical counterparts. Further, a region in which the 
VA predicts localized modes that are stable as per the 
VK criterion (hence they are stable indeed, as confirmed 
by the aforementioned numerical results) can be found 
from the variational equation dL/dW = 0: 



(1 + e-^Ag) ^ < N/{2n) < 1 



(5) 



{N = 2tt is the prediction of the VA for the norm of the 
Townes soliton) [28]). Note that the smallest value of N 
in Eq. ^ corresponds to width W = 2, which coincides 
with the radius of the circle, the stable branches of the 
p.{N) curves in Fig. [ijc) corresponding to W < 2, i.e., 
the stable modes are trapped inside the circle. Further, 
it is seen from Eq. ([5]) that, for the vanishing contrast 
between the circle and background, Ag 0, the width 
of the stability interval shrinks as AA^ « 2Tre~'^Ag « 
0.85Ag. The analysis of numerical solutions yields a very 
close result for AN. Thus we conclude that there exists 
the stable subfamily of the 2D solitons in the circle, as 
long as the SF profile is not completely flat. 

Results for the annuli are shown in Fig. 2. The funda- 
mental soliton still features a maximum at r = in Fig. 
[2l[a), in spite of the relatively large radius and depth of 
the inner hole in that case, p = 1 and Ag — 0.5. For a 
small inner radius p, the VA yields solutions in the range 
of (1 + e-^Agy' < N/{2tt) < (l - (1 + e-^)Agy\ cf. 
Eq. ([5]). However, simulations demonstrate that the fun- 
damental modes trapped in the annulus are stable only 
in a narrow interval near the left edge of this region — for 
instance, at 5.78 < N < 5.92 for p = 0.2, Ag = 0.5. At 
larger N, the modes are destabilized by perturbations 
which break the axial symmetry, driving the soliton's 
peak off the center and, eventually, causing the collapse. 
This instability (similar to that found for the annulus em- 
bedded into the linear medium (22j ) is not comprised by 
the VK criterion, which pertains solely to axisymmetric 
perturbations. 

We have also checked the dynamics of modes with em- 
bedded vorticity in this 2D setting, concluding that, as 
well as in other cases when the transverse or longitudinal 
modulation of the SF cubic nonlinearity is employed 22] , 
[H]) [11]- 01) vortex solitons cannot be stabilized against 
splitting by azimuthal perturbations (on the other hand, 
solitary vortices supported by the transverse modulation 
of the self-defocusing nonlinearity growing at r — >■ 00 can 
be easily made stable [25^]. 
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FIG. 2: (a) The profile of a fundamental soliton trapped in 
the annulus with p = 1, Ag = 0.5 and n = -0.216, N = 9.82. 
(b) The propagation constant of numerically found families of 
the trapped modes vs. A'^ in annuli with nonlinearity contrast 
Ag — 0.5 and different inner radii p. (c) The same as in (b), 
but as produced by the VA. 

The ID counterpart of the critical collapse occurs 
in the NLS equation with the quintic term, i(j)z = 
-{l/2)(j)^^ - [H m, [H. In optics, this model 

can be realized experimentally in colloids [3^. The quintic 
equation gives rise to the ID version of the Townes soli- 
ton, <j) = (— 3^)"^'''* e^*''^y^sech~(y^5^^, with the con- 
stant total power: A^o = \/3/87r for any /i < 0. The sim- 
plest ID counterpart of the 2D setting considered above 



is provided by the spatially modulated coefficient in front 
of the quintic term: g(x) = I -|- Ag ■ S{x), with Ag > 
[for a similar radial modulation of the cubic SF nonlin- 
earity in the 2D model, g{r) — go5 {r ~ tq), exact solu- 
tions for trapped modes can be found, but they are un- 
stable against azimuthal perturbations [131] • Then, the 
ID Townes solitons are replaced by exact solutions (f) — 

{-i^Jif'^ e-v^^sech(y=^(|x| -f 0)": with ^ defined by 

relation sinh (2-y/— = ZAg^J—2^, and total power 

iViD = V6 tan-i I \j y'-lS {Agf ^ + 1 + 3AgV^| 

taking values ■y/3/87r < A^id < -\/3/27r. The entire fam- 
ily satisfies the VK stability condition, dN/dfi < 0, and 
is indeed stable. 

In conclusion, it is demonstrated that the axisymmet- 
ric modulation of the strength of the SF cubic nonlin- 
earity with sharp edges stabilizes a family of 2D solitons 
against the critical collapse. For small contrast Ag, the 
width of the stability interval of the soliton's norm scales 
as AN ^ Ag. For the annulus, the stability interval is 
strongly reduced by azimuthal perturbations. 
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